Non-local conservation in the coupling field: effect on critical dynamics 
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We consider the critical dynamics of a system with an n-component non-conserved order pa- 
rameter coupled to a conserved field with long range diffusion. An exponent a characterizes the 
long range transport, a = 2 being the known locally conserved case. With renormalisation group 
calculations done upto one loop order, several regions are found with different values of the dynamic 
exponent z in the a — n plane. For n < 4, there are three regimes, I: nonuniversal, a dependent z, 
II: universal with z depending on n and III': conservation law irrelevant, z being equal to that in 
the nonconserved case. The known locally conserved case belongs to regions I and II. 
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Conservation laws play significant roles in critical dy- 
namics in thermodynamic systems. The dynamics of 
such systems have been studied over the last several 
decades 0-^]. These conservations, as constraints, gener- 
ally distinguish the various dynamic universality classes 
for a given equilibrium class. For example, in Ising fer- 
romagnets the order parameter may not be conserved 
(model A of ref. ||) though it is so for a binary mixture 
(model B). In antiferromagnets or binary alloys under- 
going an order-disorder transformation, a secondary pa- 
rameter, the total magnetization, is conserved (model C) 
. These classes are characterized by the dynamic ex- 
ponent z that connects the divergence of the time (r) and 
length (£) scales near the critical point, r ~ £ z . For sev- 
eral cases renormalisation group arguments give 
exact result for the dynamic exponent z. However, the 
numerical studies give somewhat different values which 
do not always agree with each other ||fj] . 

Most of the earlier studies assume local conservation 
of the conserved field, though local conservation is not a 
necessary, but a sufficient, condition for conservation in 
general. Recently the issue of non-local conservation has 
been addressed in the case of conserved order parameter 
(model B) either by relaxing local diffusion or allowing in- 
finite range transport |^|-[l3[ ■ For the nonconserved case, 
numerical simulations have shown that the dynamic ex- 
ponent can depend on the nature of the dynamics also 
jl4]]. It is therefore important to determine the robust- 
ness of the wellknown dynamic universality classes. 

Allowing long range diffusion in the Ising model means 
that the Kawasaki spin exchange now takes place with 
further neighbour spins also. The case where the range 
of the exchange becomes infinite corresponds to what we 
shall call "global" conservation. For the conserved order 
parameter case, model B, the effect of global conserva- 
tion seems to be drastic. In fact the dynamic correlation 
functions of the q state Potts model has been shown to be 



analogous to that the Ising model with fixed magnetisa- 
tion and long range interactions (i.e., essentially model B 
with non-local conservation) ||To|| . Dynamics for local and 
global conservation are also expected to be different [^3) . 
For O(n) symmetric order parameter, the conservation 
law has been found to be irrelevant for global conserved 
case (i.e. infinite range Kawasaki dynamics) and 
the dynamic exponent z is argued to be the same as that 
of model A. The situation is expected to be more subtle 
when a conserved quantity is coupled to a nonconserved 
order parameter as for example in model C. 

In the case of model C, the dynamics of the non- 
conserved order parameter is slowed down compared to 
model A because of its coupling to a scalar nonordering 
conserved field. Even though the conserved field is not 
critical, its coupling to the order parameter induces long 
range correlations. These singular correlations in turn 
affect the relaxational dynamics, yielding za < zc < zb 
where the subscript denotes the model name. With local 
conservation, model C has several regions in the d — n 
plane (where d is the spatial dimension and n the num- 
ber of components of the order parameter). There are 
three different regions corresponding to different values 
of the exponent z. Above n = 4, the behaviour is model 
A-like with z = 2 + crj, i.e, there is correction to 0(e 2 ), 1 
where e = 4 — d . For 2 < n < 4, the values of z are 
2+2a/(nv), and z = 2-\-a/v for n < 2, where a and v are 
respectively the specific heat and the correlation length 
exponents. (These values of z are for a > 0, otherwise 
one takes a — and z once again becomes za) The early 
time behaviour of dynamical models with short range ini- 
tial correlations also gives a new exponent p5[-p^|. Now, 
the emergence of the long range correlation in the con- 
served field, the hallmark of model C, is an equilibrium 
property. Then, are the universal dynamic behaviours, 
as reflected in the d — n plane, sensitive to the nature 
of the dynamics? We study this problem by considering 



x The equilibrium exponents a, v and r\ are obtained from 
the scaling of the specific heat (~ |t"|~ q ), correlation length 
(~ \r\~"), and correlation function (~ |x|~ (d ~ 2+,7) at r = 0) 
respectively, where r is the deviation from a critical point. 
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a long range transport for the conserved field, a conve- 
nient way of interpolation between the local to the global 
conservation. 

Let us consider the dynamics of an n component order 
parameter field <f> coupled to a scalar noncritical but con- 
served field m which may be the energy density or the 
magnetisation in the Ising antiferromagnetic case or par- 
ticle density in the binary alloy case, or annealed mobile 
impurities. We allow for subdiffusive or superdiffusive 
transport of the conserved field. This means long range 
transport for particles or the impurities Jl9| or long range 
exchange of spins in magnetic systems. The Hamiltonian 
for model C reads, 

H = J d d x[ir0 2 (x) + i(V0(x)) 2 + ^ 4 (x) 

+ 7 ^(xWx) + ic-V(x)], (1) 

where, 2 (x) = £^ 4 (x) = EIU tf(x)]V = 

is the mean field critical point, and C > ensures 
noncriticality of m. 

For purely relaxational dynamics, the equation of mo- 
tion obeyed by <f> is 



1 d(j) a (x,t) 

r ot 



OH 



(x,i) 



+ f]a(^,t), (2a) 



while to satisfies the following dynamical equation in 
Fourier space 



P(k) 



Om(k, t) 



OH 



Ot dm(-k,t) 
with the noise obeying 



■C(M) (2b) 



(0=0, ( Va ) = 0, (3a) 
<t? (x, t)r] a , (x', 0) = 2r<5(x - x')S(t - i')<W (3b) 
(C(k, t)C(k', 0) = 2P(k)S(k + k')S(t - t'). (3c) 

For model C with local conservation one chooses, P(k) = 
l/(Xk 2 ). Here we use the general form H 



P(k) 



1 



l Kk° 



(3d) 



For conservation to hold one strictly needs to satisfy 
the continuity equation only. For local conservation, 
the additional assumption is that the current is deter- 
mined by the local change in the chemical potential (this 
is the case when a = 2 and Fick's law is obeyed, i.e, 
j(x) oc V/u(x), where j is the current density and /j(x) is 
the chemical potential in the particle picture). When the 
conservation is non-local, Fick's law is no longer valid. 
This is analogous to the invalidity of Ohm's law for dif- 
fusing electrons following Levy flight instead of nearest 
neighbour hopping pQ| . In the long range diffusion case, 
instead of Fick's law, one can have a more general consti- 
tutive equation j(x) oc J /c(x— x')V/i(x')dx' (which gives 



back Fick's law for k — <5(x — x')). This for a power law 
decay of k(x) on Fourier transformation gives Eq. ( pb| ) 
with Eq. (3d). We further include the momentum inde- 
pendent term, a kinetic coefficient, Ao which is generated 
by the RG itself. A value of a 7^ 2 simulates the non- 
locality. The only constraint on a to ensure conservation 
is a > 0. The limit a — > + is the case corresponding to 
global conservation. 

We use the momentum shell renormalisation group 
method where the Fourier components 4>k (t) for the hard 
modes A/6 < k < A are eliminated and the scale changes 
k = k/b and t = b z t are made. We eventually set A = 1. 
The long distance long time results are obtained from the 
flow of the renormalized parameters. The static fixed 
points (at one loop) of this model are easily obtained 
(e = 4 - d) H as follows: 



* In + 2 

r = e 

2n + 8 



(here u = u 



AK d {n + 8) 
i 7 2 C, and K d = 2 1 - d 7r- d / 2 r(d/2)). 



d 2 cy = 



2ni>Kd 



And either 7* = or 



afv = (4 - n)e/(n + : 



(4a) 
(4b) 

(4c) 
(4d) 



a is the exponent of the temperature dependence of the 
largest term in the specific heat, i.e, a — a if a > 
and a = otherwise (a is the specific heat exponent). 
For 71 > 4, 7* = is the only fixed point. In general, 
7* = is the fixed point for n > n c where n c also has an 
e expansion. 

Considering the various vertex functions, the renor- 
malisation equation for the transport and kinetic coeffi- 
cients can be obtained || . The RG transformation being 
analytic cannot renormalise the singular terms of P(k), 
as in the local case. However the transport coefficient 
generates (and renormalizes) the kinetic coefficient, A , 
and also renormalizes the kinetic coefficient of the or- 
der parameter. Defining the dimcnsionless parameter 
/Jer = rC/Acr, the recursion relations at the static fixed 
point of Eq. (fk]) can be written as (choosing thin shells 
b = 1 + 51 with SI -> 0) 



dT- 



01 



= r- 1 \-z + 2- 



2a /vn 



9A" 1 



01 
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(5a) 

(5b) 
(5c) 



Here \x = u a + /io, with u = K^-. The above one loop 
contributions in the dynamic quantities come from the 
7 vertex. (It maybe noted that we have excluded 77 as 
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7/ = to first order in e. This is not crucial as there is a 
cancellation in model C in higher order terms Js|,§].) 

The dynamic exponent z is chosen to keep T scale in- 
variant. Eq. (pal) then gives 



z = 2 



2a /vn 



(6) 



+ 

The fixed point of is obtained from eqs. (^) and ([Ic]) 

a/v 



Mo 



2(z - a/i/) 



The fixed points of fi a are therefore 0, oo or a non-zero 
finite value given by 



a - 2 + <S/V 

Z — (7 H 

1 V 71 



(J) 



Stability of fixed points and exponents: 

Let us analyse the stability of the fixed points for the 
different values of a. 

Above n = 4, 7* = is the only fixed point and there- 
fore the conserved field gets decoupled in long lengthscale 
limit. This is nothing but the model A fixed point where 
the ordering and non-ordering fields are no longer cou- 
pled and therefore z = za- To recover za = 2 + crj one 
needs to go to higher loops, where the u vertex starts 
contributing. 

= is the fixed point valid for a < 2 — a/v and the 
exponent z = 2 + 0(e 2 ) for all n . At this fixed point the 
conservation of the coupling field simply becomes irrel- 
evant. One therefore again recovers model A exponent 
z = za, although 7* ^ 0. 

We note that as /i* cannot be negative and a ~ O(e), 
(fjj) will give a physical fixed point only if n < 2/(1 +p) 
where p = (cr — 2) /(a/v) and p > — 1. For 4 > n > 
2/(1 +p), p > — 1, Ha = °o will be the stable fixed point. 
The latter gives n = 4 at p = —0.5, so for —1 < p < —0.5, 
we find that a non-zero finite fixed point for /1* is always 
stable. In particular, for a — 2, we get back the results of 
||: z = 2 + a/v for n < 2, z = 2 + 2a/ nv for 2 < n < 4 
and z = 2 + 0(e 2 ) for n > 4. 

Thus we have the following regions in the a — n plane: 

I. z — cr + a/v. 
For 

a = 2+pa/v,n < 2/(1 + p); p> -1 

II. z = 2 + 2a/nv. 
For 

(7 = 2+pa/v,4 > n > 2/(1 + p); p > -1/2 

III. z = 2 + 0(e 2 ), 7* = 0. 
For all a, n > 4. 

Ill' z = 2 + 0(e 2 ), 7* ^ 0. 



For a < 2 — a/v (i.e., p — — 1), n < 4. 

These are shown schematically in Fig. 1. We may add 
that there are discontinuities of z only at the boundaries 
(7 = 2 — a/v and n = 4. The model C like regions for 
(7 = 2 are in regions I and II with nothing special for 
the cr = 2 line. The value of z is also continuous across 
the boundary of regions I and II. Of the two regions, I 
is nonuniversal with a a dependent dynamic exponent 
z = cr + a/v while II is robust and independent of a. In 
region II, z is dependent on n explicitly while in region I, 
n dependence is through the equilibrium exponents only. 
As one crosses the a = 2 line, region II extends towards 
smaller values of n for a > 2 and the opposite happens 
for a < 2 where it vanishes at n = 4, a = 2 — a/2v and 
region I survives as long as a > 2 — a/v. 

The general a case, in addition, allows for a new 
regime where the conservation law becomes irrelevant 
even though the coupling survives (7 7^ 0) in the long 
length scale limit and induces long range correlations in 
equilibrium. The global conservation case belongs to this 
class. 

Since the nonrenormalizability of \& is valid to all or- 
ders, we expect, as in the local model C, the dynamic 
exponent z, as obtained here, for regions I and II to be 
exact. The above analysis is best for d close to 4 as we 
have done calculations upto O(e) only. For d close to 4, 
in fact, region I exists for values of cr > 2 and a < 2 and 
very close to 2 only as a ~ O(e). Continuing with this 
result in lower dimensions, one expects region I to exist 
over a larger range of a(> 2) so that the line separating 
regions I and II in the positive half of the a — 2 axis in 
Fig. 1 goes higher up for lower dimensions. However, the 
line separating regions II and III (or III' and III) in Fig 1 
is shifted towards lower values of n for lower dimensions 
as a becomes negative at a critical value of n — n c where 
n c also has an e expansion. (Fig. 1 has been drawn for 
d close to 4). In the d — n plane therefore, region I is 
in general depleted for a value of a higher than 2. This 
depletion is maximum for d — 4. For a < 2 — a/v, only 
region III or III' survive with z = za for all d. 

Some of our results are in contrast with the long range 
diffusion in model B, where z is explicitly dependent on 
a as long as the conservation is relevant. The dynamical 
exponent is in general expected to increase for a system 
with constraint, and therefore z > za- However, global 
conservation (i.e., a — > + ) becomes irrelevant and one 
recovers z = za in this limit j|,[l2|,[ll| . This would the 
case for a < a c = za — 2 + rj = 0(e 2 ) for model B. 
In model C, we find region III', where conservation is 
irrelevant, to be significantly larger with a c = 2 — a/v. 
Even when the conservation is relevant, we find a strongly 
universal regime, region II, where the conservation slows 
down the dynamics but the nature of the dynamics is 
irrelevant (cr independence). 

It is possible to implement this long range transport 
in a simulation. The infinite range would correspond 
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to (7 — > while any other a can be implemented by 
choosing the two spins to be flipped, in the Kawasaki dy- 
namics, at a distance I with probability P(l) ~ l~( r+cr ) 
for large I. However, it is fairly easy to identify the irrel- 
evance of the conservation in model B numerically when 
the model A exponent is recovered, as za and zb are eas- 
ily distinguishable even with moderate level of precision 
in calculations. For model C, very high precision simula- 
tions are required to verify our RG based scenario as the 
model C exponents in region I and II are different from 
that of model A by a much lesser margin. 

To summarise our results, we find the effect of long 
range diffusion in model C with an n-component or- 
der parameter to be highly significant. Characterising 
the long range diffusion by an exponent er > 0, sev- 
eral regimes with different dynamical exponents z are 
obtained in the er — n plane. For n < 4, a region with a 
a dependent z is found to exist as well as a region with 
strong universality where z is independent of a. Most 
remarkable is the existence of a region for a greater than 
a critical value where the coupling with the field survives 
in the large length scale limit but the conservation be- 
comes irrelevant and z — za ■ Such a region is also present 
in model B with non-locally conserved order parameter, 
although for a much smaller range of a. For n > 4, 
the coupling with the field is irrelevant and z — za as 
expected. The known results for the locally conserved 
model C are recovered with a = 2. 
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FIG. 1. Regions corresponding to different values of z are 
shown for p = and n. The flow diagrams in fx a — fio plane 
are also shown. The flow is towards a finite non-zero value of 



ficr in region I, towards \i a 
region IIP. 



oo in region II and \i a — in 
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